The nuclear reaction waiting points, Mg, 6 Si, S, and 3 Ar, and 
bolometrically double peaked type I X-ray bursts 
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ABSTRACT 

Type I X-ray bursts with a double peak in the bolometric luminosity have been observed from 
several sources. The separation between the two peaks are on the order of a few seconds. We 
propose a nuclear waiting point impedance in the thermonuclear reaction flow to explain these 
observations. Nuclear structure information suggests the potential waiting points: 26 Mg, 26 Si, 
30 S, and 34 Ar, which arise in conditions, where a further reaction flow has to await a /3 + -decay, 
because the (a,p)-reaction is too weak to overcome the target Coulomb-barrier and the (p, 7)- 
reaction is quenched by photo-disintegration at the burst temperature. The conclusion is that the 
effects of the experimentally unknown 30 S(a,p) 33 Cl and 34 Ar(a,p) 37 K might be directly visible 
in the observation of X-ray burst light curves. 

Subject headings: X-rays: bursts — nuclear reactions, nucleosynthesis, abundances — stars: neutron 



1. Introduction 

A type I X-ray burst results when a thermal in- 
stability in the electron degenerate atmosphere of 
a neutron star causes a thermonuclear runaway in 
the accreted matter. The sudden release of nuclear 
energy heats the atmosphere rapidly, increasing 
the luminosity a few orders of magnitude above 
the persistent luminosity within a few seconds, af- 
ter which the luminosity decreases exponentially 
as the atmosphere cools again. While most of 
these bursts show a single peak in the luminos- 
ity, some burst observations have shown a dou- 
ble peaked structure in the bolometric light curve 
(Sztajno et al. (1985); van Paradijs et al. (1986); 
Penninx et al. (1989); Seon & Min (2000) and 
Smale (2001)). Several theories have been pro- 
posed to explain such double peaks: Heat trans- 
port impedance in a two-zone model (Regev & 
Livio 1984), interactions with the accretion disc 
(Melia & Zylstra (1992)), a result of a H flash de- 
veloping into a combined H/He flash (Hanawa & 



Fujimoto 1984), shear instabilities (Fujimoto et al. 
1988) or increased proton capture on heavier nu- 
clei (Ayasli & Joss 1982). 

In this letter we propose that an interaction be- 
tween the helium flash shell and a waiting point 
impedance in the rp-process of the shells above 
it explains burst shape structure spanning over a 
time-scale of a few seconds. Since a self-consistent 
one-dimensional model attains a limit cycle equi- 
librium with repeated self-similar bursts, compar- 
isson with observations will also place restrictions 
on the accretion rate. 

2. The first burst peak: A helium flash 

Observationally burst light curves differ greatly 
in shape but their understanding are theoretically 
divided into three cases, depending on whether the 
ignition and subsequent runaway occurs in a hy- 
drogen rich shell, a helium rich shell, or in a mixed 
H/Hc shell (Fujimoto et al. 1981). Burst sources 
are situated in low mass binary systems and gen- 
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crally accrete a mixture of hydrogen and helium. 
Since the atmosphere of the neutron star has a 
temperature of Tg ~ 0.1 — 0.2, the /3-limited hot- 
CNO cycle (Wallace & Woosley 1981) converts the 
hydrogen into helium. 

For a constant accreting rate 1 of M — 0.045Me<;. 
it requires about half a day before the continuous 
accretion process has pressed the matter down 
to a depth, where the helium which has now been 
formed in significant quantities becomes thermally 
unstable. The accumulated helium at the base 
burns completely within tenths of a second and 
establishes a temperature gradient which is suffi- 
ciently large to allow a highly efficient convective 
heat transport to the surface. Observationally 
this results in a rise time of the surface luminosity 
which is less than a second (Wallace ct al. 1982) 
and it is the cause of the first peak. 

3. The second burst peak: A nuclear wait- 
ing point impedance 

The helium burning region radiates and con- 
vects heat upwards so that hydrogen which has 
not yet had the time to be fully converted into he- 
lium ignites as well and burns via the rp-process 
(Wallace & Woosley 1981; van Wormcr ct al. 1994; 
Schatz et al. 1998). Following the breakout of 
the hot-CNO cycle the reaction flow in the sd- 
shell up to 40 Ca is determined by a competition 
between the (a,p) process (Wallace & Woosley 
1981), which is very temperature dependent, be- 
cause of the strong Coloumb-barrier against alpha- 
capture, and the rp-process. The latter is only 
temperature dependent in the sense that the tem- 
perature determines the ratio of the proton-rich 
isotopes within an isotone chain and thus which 
isotopes /3 + -decay and move the flow between the 
isotone chains (Rembges et al. 1997). It is the 
detailed behavior of this flow which determines 
the burst behavior following the first peak. If the 
reaction flow is stopped at a "waiting-point" iso- 
tope, the nuclear energy release rate decreases and 
causes a dip in the surface luminosity. The surface 
luminosity rises again, when the nuclear energy re- 
lease rate - following the waiting - increases again. 

The Eddington critical accretion rate equals 1.1 ■ 10 18 g/s 
(1.7 • lO _8 M0/yr) and assumes temperature independent 
Thompson scattering in an atmosphere with a solar com- 
position. 



The /3 + -unstable even-even T z = — 1 isotopes, 
22 Mg, 26 Si, 30 S, and 34 Ar, are potential waiting 
points, which are important, because most of the 
reaction flow passes through them. The reason 
is that increasing temperatures shift the flow to- 
wards the proton driplinc, whence only two possi- 
ble flow paths circumventing them need considera- 
tion (see fig. 1). The first path depends on branch- 
ing between the /3 + -decay (Ti/ 2 ~ few seconds) of 
the even-A-odd-Z T z = (N - Z)/2 = -1/2 iso- 
tope and the proton-capture going to the poten- 
tial waiting point isotope, but because of the high 
temperature the proton-capture rate is very much 
larger, which makes the decay insignificant in com- 
parisson. The second path would depend on dou- 
ble /3 + -decays of very proton-rich isotopes between 
successive odd-Z isotones. These decays proceed 
through the even-Z isotones where the (p, 7)-(7, p) 
reactions generally shift the equilibrium towards 
the waiting points. 

Since 22 Mg, 30 S, and 34 Ar have low (p,j) Q- 
values (colored red on fig. 1) they will enter (p, 7)- 
(7,p) equilibrium with 23 Al, 27 P, 31 S, and 35 K 
respectively (Thielemann ct al. 1994), whence 
the ratio of the target and the compound iso- 
tope is only determined by a Boltzmann factor, 
exp(— ksT/Q), along with a spin statistics factor. 
For 22 Mg (Q = 126 keV), 30 S (Q = 296 keV), 
and 34 Ar (Q = 78 keV) the Q-values are suf- 
ficiently low to shift the ratio towards the for 
photo-disintegration to prevent a significant flow 
through the proton-capture compound isotope. 
This is not the case for 26 Si (Q = 859kcV) where 
the equilibrium balance is shifted towards 27 P, so 
unless the Q- value has been overestimated, 26 Si 
is not a waiting point, because it is destroyed via 
the /3+-decay of 27 P. 

In addition to the (p, "/)-("/, p) equilibrium 
and the /3 + -decay a third branching of the flow 
from these isotopes originates in the (a,p)- 
process which commences on 14 O and proceeds 
via 14 0(a,p)(p, 7 ) 18 Nc(a,p)(p, 7 ) 22 Mg(a,p)(p, 7) 
26 Si(a,p)(p, 7 ) 30 S(a,p)(p, 7 ) 34 Ar(a,p)(p, 7) 38 K. 

An isotope is a potential waiting points if the 
(a,p)-reaction is so weak that the major part of 
flow must await /3 + -decay. Fig. 2 shows a plot of 
the condition 

X 0+ = N A Y a A(T) < av > {a , p) , (1) 

which demarcates the (p, T)-space between a dom- 
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Fig. 1. — Targets with (p, 7) Q-values of < 0.0 
(white), 0.0-0.5 (red), 0.5-2.5 (yellow), 2.5-3.5 
(green), 3.5-5.5 (blue), 5.5-8.0 and (purple). Sta- 
ble nuclei are color-coded grey. The reaction flow 
rates, fa = Yi^j — Y,-_>i, during the burst peak 
temperature are shown in solid lines. The thick- 
ness indicates the strength of the reaction flow 
with the exception of isotopes in (p, 7)-(7,_p) equi- 
librium which are shown in thick lines though the 
net flow is close to zero. The reaction paths which 
circumvent the waiting points are indicated with 
dashed arrows. 

inant /3 + -decay (below the line) and a dominant 
proton capture (above the line) assuming that 
Y a = 0.5/4 = 0.125. The plot ignores the (p,j)- 
reaction leak which would move the line down- 
wards. The temperature dependent function A(T) 
accounts for the uncertainty in the reaction rate 
which in this paper is set to be constant for the 
sake of simplicity. Plotting eq. [1] with A(T) = 
10 2 and A(T) = 10~ 2 gives an estimate of the 
effect of the uncertainty of the reaction rate. 

From fig. 2 it is seen that 22 Mg(a,p) 25 Al dom- 
inates 22 Mg(/3+) 22 Na even if 22 Mg(a,p) 25 Al had 
been underestimated by a factor 100. It is there- 
fore unlikely that 22 Mg is a waiting point. How- 
ever, both 30 S(a,p) 33 Cl and 34 Ar(a,p) 37 K block 
the flow during the burst conditions, because the 
temperature never rises sufficiently to break the 
Coulomb barrier. 

Therefore the destruction rate (see e.g. van 
Wormer et al. (1994)) of the waiting point iso- 
topes depends on the interplay between three rates 
whose exact details depend on nuclear data and 
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Fig. 2. — The demarcation lines of eq. [1] are 
plotted thrice with A(T) = {10" 2 ,1,10 2 } and 
Y a = 0.125 for the 22 Mg, 26 Si, 30 S, and 34 Ar 
isotopes. The demarcation line for underestimat- 
ing 34 Ar(a,p) 37 K by a factor 10 2 falls outside the 
graph. In addition the solid line shows the thermo- 
dynamic trace of the shell conditions during a full 
burst cycle which is computed in section 4. The 
arrows indicate the time direction of the cycle. 

the density and temperature through the interac- 
tion with the hydrodynamics. 

4. Computational results and conclusion 

The combined effect of the two peaks has 
been computed using a new self-consistent one- 
dimensional model. It is based on modified ver- 
sion of AGILE (Liebendorfer et al. 2002) with 
added turbulent mixing using a diffusion equation 
and convective heat transport based on Mixing 
Length Theory. The code tracks energy trans- 
port to a high precision and realistically considers 
the heat transport between the atmosphere and 
the core. The reaction network, which includes 
proton-, neutron-, and alpha-induced reactions 
and accounts for neutrino losses, comprises 298 
isotopes extending up to Tellurium. Further tech- 
nical details will be described in an upcoming 
paper (Fisker et al. in preparation). 

Having eliminated 22 Mg and 26 Si, we inves- 
tigate 30 S and 34 Ar by assuming a local rest 
mass accretion rate of M = 5 • 10 16 g/s and 
computing the bolometric luminosity curves for 
the current reaction-rates which are based on the 
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code NON-SMOKER (Rauscher & Thielemann 2001). 
In addition luminosity curves were computed as- 
suming that the 30 S(a,p) 33 Cl and 34 Ar(a,p) 37 K- 
reaction rates have been underestimated by a fac- 
tor 100 i.e. A(T) = 100. This implicitly tests any 
uncertainty in the /?+- decay half lives. The lu- 
minosity curves are shown in fig. 3, while a trace 
of the temperature and density is shown in fig. 2 
using the unaltered reaction rates. Fig. 3 shows 
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Fig. 3. — The computed luminosity is plotted as 
a function of time for the three cases. The three 
luminosity curves have been synchronized to make 
the second burst peaks coincide. 

that if the 30 S(a,p) 33 Cl-reaction were faster then 
the dip in the luminosity would decreased or pos- 
sibly be eliminated. Increasing the 34 Ar(a,p) 37 K- 
reaction decreases the dip slightly indicating that 
the 34 Ar-waiting point plays a smaller role in cre- 
ating the full dip in the unaltered model. In both 
cases increasing the rate of either (a,p)-reaction 
increases the luminosity of the second peak indi- 
cating that a combination would increase further. 

Convective heat transport defines the rise of the 
first burst peak, whence a time dependent model 
of the turbulent part of the kinetic energy could 
change the shape or magnitude of the first burst 
peak. In addition a more effective convective heat 
would increase the temperature in the upper re- 
gions which would make the (a,p)-reaction flow 
stronger despite the resulting drop in density from 
the temperature increase and vice versa. 

In conclusion the most uncertain factors for 
the impedance driven double peak is the model 
of convection, the accretion rate, the restriction 



to a spherically symmetric model and the uncer- 
tainty in the theoretical (a,p) rates 30 S and 34 Ar. 
However, even if the effect may not be decidedly 
observable in the light curve, the waiting point 
impedance of 30 S and 34 Ar and possibly 22 Mg and 
26 Si is inherent to the system and an experimental 
confirmation of the theoretical rates would help 
to determine the limits of the fast flow through 
the sd-shell isotopes, whereas improved hydrody- 
namics is important to affirm the (p, T)-parameter 
space of the (a,p)-reactions. 

JLF and FKT acknowledge support from the 
Swiss NSF grant 2000-061822. 
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